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ABSTRACT 

We  describe  multi-baseline  observations  of  a  geostationary  satellite  using  the  Navy  Precision  Optical  Inter¬ 
ferometer  (NPOI)  during  the  glint  season  of  March  2015.  We  succeeded  in  detecting  DirecTV-7S  with  an 
interferometer  baseline  length  of  8.8  m  on  two  nights,  with  a  brief  simultaneous  detection  at  9.8  m  baseline 
length  on  the  second  night.  These  baseline  lengths  correspond  to  a  resolution  of  ~4  m  at  geostationary  altitude. 
This  is  the  hrst  multiple-baseline  interferometric  detection  of  a  satellite. 
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1.  INTRODUCTION 

Developing  the  ability  to  resolve  geostationary  satellites  (geosats)  from  the  ground  for  Space  Situational  Aware¬ 
ness  is  currently  drawing  wide  interest.  The  largest  single  telescopes  are  just  capable  of  resolving  the  large-scale 
structure  of  these  satellites,^  but  resolving  details  requires  the  use  of  other  techniques,  such  as  optical  inter¬ 
ferometry. 

The  Navy  Precision  Optical  Interferometer  (NPOI)^’^  is  the  hrst,  and  so  far  only,  interferometer  to  succeed 
in  detecting  a  satellite.  We  observed  DirecTV-9S  with  a  16  m  baseline  during  the  glint  season  of  March 
2008,  with  followup  observations  in  2009'^  resulting  in  a  glint  model  consisting  of  a  resolved  (>  3.7  m)  and  an 
unresolved  (~  1  m)  component,  with  roughly  equal  huxes. 

One  of  the  lessons  from  those  observations  is  that  the  baseline  length  B  of  16  m,  the  shortest  available 
at  the  NPOI  at  the  time,  is  somewhat  too  long  for  observing  these  targets,  for  which  the  typical  size  scale  is 
several  meters.  The  A  =  556  nm  to  845  nm  wavelength  range  used  in  those  observations  produced  a  resolution 
0res  «  A/i?  «  35  nrad  to  50  nrad,  corresponding  to  1.3  to  2  m  at  geostationary  altitude.  Structure  at  that 
scale  produces  a  fringe  contrast  V  ~  0.2,  and  structure  that  is  only  20%  larger  than  6,-es  produces  no  fringe 
contrast  at  all.  At  still  larger  scales,  the  fringe  contrast  recovers,  but  only  weakly. 

The  lack  of  shorter  baselines  puts  us  at  a  disadvantage  in  two  ways.  First,  it  makes  it  difficult  or  impossible 
to  detect  and  track  fringes  in  real  time  and  thereby  adjust  the  internal  optical  paths  in  the  interferometer, 
particularly  because  atmospheric  turbulence  over  the  array  elements  forces  us  to  detect  fringes  and  readjust 
optical  paths  on  20  ms  timescales.  Second,  if  we  have  data  only  from  long  baselines,  those  data  do  a  poor  job  of 
characterizing  the  larger-scale  structure  of  the  target.  A  better  representation  of  the  larger  scales  representative 
of  the  satellite  bus  requires  data  from  shorter  baselines,  or  in  interferometric  parlance,  sampling  smaller  u  —  v 
spacings,  where  u  and  v  are  the  spatial  frequencies  i?Ew/A  and  Bns/A  sampled  toward  the  target  by  the 
east-west  and  north-south  components  of  the  baseline. 

Since  2009,  the  NPOI  has  brought  several  new  array  stations  into  operation  that  make  shorter  baselines 
available  (see  Fig.  1).  In  particular,  we  now  have  an  8.8  m  baseline  between  stations  W4  and  AC,  which  yields 
resolutions  of  2.3  to  3.5  m  at  geostationary  altitude,  and  a  9.8  m  baseline  between  stations  AC  and  E3.  Having 
two  short  baselines  that  share  an  array  element,  such  as  these,  makes  it  possible  to  “bootstrap”  them  together 
to  make  a  third  baseline,  W4  to  E3,  of  18.6  m.  If  fringes  are  detected  and  tracked  on  the  two  shorter  baselines. 
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Figure  1.  Layout  of  the  inner  stations  of  the  NPOI.  The  black  squares  are  the  astrometric  stations  AW,  AC,  AE  and 
AN.  The  black  triangles  show  the  new  imaging  stations,  W04,  EOS  and  NOS,  commissioned  for  the  observation  of  geosat 
glints.  The  black  lines  show'  the  baselines  W04-AC  and  AC-EOS  used  in  these  observations  of  DirecTV-7S. 

the  third  basehne  is  also  correctly  phased,  and  it  can  take  data  without  fringes  being  detectable  on  it  in  real 
time,  a  concept  due  originally  to  Roddier.^ 

Using  three  or  more  array  elements  also  allows  us  to  obtain  closure  phases,  the  sum  of  the  fringe  phases 
aroimd  the  triangle  of  basehnes.  Because  atmospheric  perturbations  of  the  basehne  phases  cancel  in  the  sum,® 
these  closure  phases  represent  some  (though  not  all)  of  the  phase  information  of  the  image  of  the  target.  In 
particular,  they  can  be  used  to  better  determine  the  relative  locations  of  the  components  of  the  satellite. 

2.  OBSERVATIONS 

2.1  NPOI 

The  12  cm  apertmes  of  the  NPOI’s  array  elements,  in  conjunction  with  atmospherically-induced  fringe  motion 
that  requires  adjusting  the  internal  path  lengths  at  ~50  Hz,  limit  the  NPOI’s  sensitivity  to  visual-band 
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magnitudes  V  ^  6  mag.  During  most  of  the  year,  most  geosats  are  much  fainter,  typically  ~  10  mag  to 
15  mag.  Twice  yearly,  however,  some  satellites  whose  solar-panel  arrays  are  oriented  north-south  can  glint 
brightly  by  specular  reflection  of  sunlight.  Because  we  see  the  Sun  reflected  by  the  panels,  the  “glint  season” 
occurs  when  the  apparent  declination  of  the  Sun  is  the  same  as  that  of  the  satellites.  From  the  NPOI,  the 
geosat  belt  is  at  an  apparent  declination  of  —6°,  so  the  glint  seasons,  which  are  about  10  days  long,  are  centered 
on  March  1  and  October  11,  when  the  Sun  is  south  of  the  Equator.  Glints  can  be  as  bright  as  1®*  magnitude, 
and  typically  last  for  5  to  20  min  each  night. 

We  observed  DirecTV-7S  with  the  NPOI  on  the  nights  of  2015  March  5  and  6,  using  stations  W4,  AC, 
E3,  and  N3.  We  used  four,  rather  than  the  maximum  of  six,  stations  in  part  because  having  fewer  baselines, 
and  hence  fewer  interference  fringes,  in  the  data  makes  detecting  the  fringes  easier.  In  the  event,  N3  had 
a  malfunction,  leaving  us  with  three  stations  and  three  baselines.  Eor  most  of  the  February  26  to  March  6 
observing  run,  the  array  configuration  was  purely  academic  due  to  heavy  snowfall,  and  although  the  skies 
cleared  for  the  5*^  and  G**',  we  had  very  poor  seeing  (high  atmospheric  turbulence). 

The  NPOI  “Classic”  beam  combiner^  used  for  these  observations  uses  pathlength  modulation  to  scan  the 
interference  fringes.  The  feed  beam  from  each  array  element  traverses  a  delay  line,  which  imposes  a  500  HZ 
triangle  wave  modulation  on  each  beam.  Modulation  amplitudes  are  multiples  of  1  /rm,  creating  triangle- 
wave  optical  path  difference  (OPD)  variations  with  amplitudes  of  1,  2,  or  3  /rm  between  beams  when  they 
are  combined.  We  used  two  of  the  Classic  combiner’s  three  output  beams,  each  fed  to  a  spectrograph  that 
records  16  spectral  channels  in  the  556-845  nm  wavelength  range.  Spectrograph  1  recorded  all  three  baselines 
simultaneously,  while  baseline  W4-AC  appeared  as  well  in  spectrograph  2.  The  Classic  combiner  collects  data 
for  an  integral  number  of  wavelengths  in  each  wavelength  channel,  so  for  OPD  amplitudes  of  1,  2,  and  3  /xm, 
the  fringe  frequencies  in  the  data  are  fc  =  1,  2,  and  3  fringes  per  modulation  period,  corresponding  to  baselines 
W4-E3,  W4-AC  and  AC-E3,  respectively. 

2.2  Supporting  observations 

Observing  geosats  with  an  interferometer  differs  from  observing  stars  in  two  ways:  we  don’t  know  precisely 
where  they  are,  and  we  don’t  have  much  time  to  find  them.  As  for  the  position,  pointing  the  individual  array 
elements  at  the  geosat  is  straightforward;  what  is  difficult  is  knowing  the  OPD  for  the  particular  target.  For 
a  9  m  baseline,  shifting  the  OPD  by  100  /xm  corresponds  to  shifting  the  sky  position  at  which  the  baseline 
is  pointed  by  2.3  arcsec,  while  the  position  calculated  from  the  satellite’s  two-line  orbital  elements  (TLEs) 
published  on  line  is  typically  10  times  less  precise.  Since  searching  100  /xm  in  OPD  space  can  take  several 
minutes,  knowing  the  satellite  position  to  1  arcsec  or  better  can  make  the  difference  between  finding  fringes  at 
the  correct  OPD  while  the  satellite  is  still  glinting  and  missing  the  glint  altogether. 

We  used  the  Naval  Observatory  Flagstaff  Station  40-inch  telescope  in  the  R  band  (A  «  660  nm)  to  measure 
the  positions  and  drift  rates  of  our  targets  until  the  glints  became  bright  enough  at  i?  «  7  to  saturate  the 
camera.  We  then  extrapolated  forward  for  the  next  hour  to  generate  positions  for  steering  the  OPDs  at  the 
NPOI.  The  positional  uncertainty  from  our  position  measurements  and  extrapolations  was  better  than  1  arcsec. 

3.  RESULTS 

3.1  Photometry 

The  16  spectral  channels  of  the  NPOI  provide  a  means  of  doing  low-precision  comparative  photometry.  The 
uncalibrated  spectrum  of  DirecTV-7S  is  presented  in  Fig.  2,  along  with  the  spectrum  of  a  G9  giant  star, 
C  Hydrae,  observed  the  same  night.  The  spectrum  of  the  star  was  divided  by  1.9  in  order  to  match  the 
spectrum  of  DirecTV-7S  at  845  nm  and  highlight  the  color  difference  between  the  star  and  the  satellite.  From 
the  V  (A550  nm)  and  I  (A805  nm)  magnitudes  of  (  Hydrae  (3.10  mag  and  1.90  mag,  respectively”^)  and  our 
measured  flux  ratios,  we  calculate  that  DirecTV-7S  magnitudes  were  V  =  4.5  mag  and  I  =  2.6  mag  during 
the  period  when  fringes  were  detected.  With  V  —  I  =  1.9  mag,  this  satellite  is  significantly  redder  than  the 
Sun  {V  —  I  =  0.70  mag®),  the  source  of  light  it  reflects.  This  result  is  consistent  with  other  observations  of 
geostationary  satellites.®^^^ 
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Figure  2.  Uncalibrated  spectrum  of  the  geosat  DirecTVTS  observed  with  NPOI  Classic  spectrograph  1  (solid  black  line) 
and  the  star  ^  Hydrae  (HR  3547;  green  dashed  line).  The  vertical  axis  shows  the  number  of  photons  detected  per  2  ms 
frame.  The  618  nm  channel  is  dead  and  has  been  set  to  zero. 
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We  obtained  fringes  on  DirecTV-7S  on  a  single  baseline,  W4-AC  on  the  night  of  2015  March  5  (UT),  but 
because  the  combiner  software  was  searching  for  a  second  baseline,  the  data  were  not  recorded.  The  next 
night,  fringe  detection  continued  to  be  intermittent  due  to  poor  weather,  but  we  were  able  to  detect  and  record 
fringes  over  a  period  of  27  s,  during  which  we  tracked  the  fringe  on  baseline  W4-AC  for  approximately  4  s  and 
the  fringe  on  baseline  AC-E3  for  a  smaller  fraction  of  the  time.  The  observations  were  further  complicated 
by  the  fact  that  they  were  done  during  the  last  days  of  the  glinting  season,  which  means  that  the  glints  were 
shorter  in  duration  and  less  extreme  in  brightness  than  they  would  have  been  at  the  peak  of  the  season. 

The  fringe  power  spectra  shown  in  Figs.  3  and  4  show  the  detection  from  March  6  in  each  of  12  spectral 
channels.  Figure  3  shows  spectrograph  2,  which  observed  only  baseline  W4-AC  at  a  fringe  scanning  frequency 
k  =  2  cycles  per  modulation  period.  The  bias-subtracted  fringe  power  (green  trace)  shows  a  strong  signal 
in  the  845  nm  channel  (chi),  with  declining  fringe  power  as  the  wavelength  gets  shorter.  Figure  4  shows 
spectrograph  1,  with  which  observed  all  three  baselines  were  observed.  Again,  there  is  a  strong  signal  at  A:  =  2 
from  baseline  W4-AC  in  channel  1,  but  there  is  also  a  weaker  signal  from  baseline  AC-F3  at  fc  =  3.  These 
fringe  power  signals  can  be  followed  to  channel  7  (700  nm)  or  8  (681  nm). 

Figure  4  does  not  show  evidence  of  a  fringe  at  A:  =  1,  corresponding  to  the  bootstrapped  W4-E3  baseline, 
due  to  some  combination  of  the  bad  weather,  the  size  of  the  target,  and  the  short  time  over  which  we  detected 
fringes,  especially  on  the  AC-E3  baseline. 


4.  FUTURE  WORK 

We  are  currently  analyzing  our  data  from  DirecTV-7S  and  from  stars  observed  during  2015  March  6  to 
determine  methods  of  eliminating  frames  taken  during  moments  when  the  fringes  were  not  present  and  of 
averaging  the  on-fringe  data.  With  appropriately  filtered  data,  we  will  be  in  a  position  to  model  the  data  and 
determine  the  sizes  and  relative  fluxes  of  the  glinting  structures. 

The  next  glint  season  at  the  NPOI  occurs  in  early  October  2015.  We  plan  to  use  the  same  compact  NPOI 
configuration.  Backend  software  and  hardware  improvements  now  under  way  may  allow  us  to  record  data  in 
32  spectral  channels  (845  nm  to  450  nm)  with  all  three  of  the  output  beams  of  the  current  Classic  combiner. 
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SPECTROGRAPH  2 


Figure  3.  Fringe  power  as  a  function  of  fringe  scanning  frequency  k  fringes  per  modulation  period  in  DirecTV-7S 
Spectrometer  2  data  from  March  6.  Each  frame  corresponds  to  a  spectral  channel,  with  wavelengths  from  845  nm 
(chi)  to  604  nm  (chl2).  This  spectrometer  observed  only  the  W4-AC  baseline,  which  had  a  fringe  scanning  frequency 
fc  =  1.  The  black  line  shows  the  power  sjjectrum  obtained  with  the  coherent  (on- fringe)  scan.  The  red  line  is  the 
power  spectrum  bias'*  obtained  using  an  incoherent  (off-fringe)  scan,  scaled  based  on  the  flux  ratios  of  the  coherent 
and  incoherent  scans.  The  green  line  is  the  bias-corrected  power  spectrum. 
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Figure  4.  Fringe  power  as  a  function  of  fringe  scanning  frequency  k  fringes  per  modulation  period  in  Spectrometer  1 
data  from  March  6.  This  spectrometer  observed  all  three  baselines,  W4-E3,  W4-AC  and  AC-E3,  with  fringe  frequencies 
A:  =  1,  2  and  3,  respectively.  As  in  3,  the  black,  red,  and  green  lines  show  coherent  (on-fringe),  bias  (off-fringe),  and 
bias-corrected  power  spectra,  respectively.  Channel  11  is  a  dead  channel. 
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